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I. INTRODUCTION 



The established evidence that neutrinos oscillate and possess a small mass below the 
electron volt (eV) scale [l[ necessitates physics beyond the Standard Model (SM), which 
could manifest itself at the CERN Large Hadron Collider (LHC) and the Fermilab Teva- 
tron, and/or in low energy experiments which search for lepton flavour violation (LFV) (2)]. 
Consequently, models of neutrino mass generation which can be probed at present and 
forthcoming experiments are of great phenomeno logical interest. 

Neutrinos may obtain mass via the vacuum expectation value (vev) of a neutral Higgs 
boson in an isospin triplet representation j3|-0]- A particularly simple implementation of this 
mechanism of neutrino mass generation is the "Higgs Triplet Model" (HTM) in which the 
SM Lagrangian is augmented solely by a SU (2) triplet of scalar particles A with hypercharge 
In the HTM, neutrinos acquire a Majorana mass given by the product of 
a triplet Yukawa coupling (hij, with i,j 



there is a direct connection between hi 



= e,/i,r) and a triplet vev (v&). Consequently, 
<,jj and the neutrino mass matrix, which gives rise 
to phenomenological predictions for processes which depend on h^. A distinctive signal 
of the HTM would be the observation of doubly charged Higgs bosons (if ±:t ) whose mass 
(Mff±±) may be of the order of the electroweak scale. Such particles could be produced 
with sizeable rates at hadron colliders via the processes qq — > H ++ H jsl— 12 1 and q'q — > 
H ±± H T ' 



13, 14 



The first searches for H ±i: at a hadron collider were carried out at the 
Fermilab Tevatron, assuming the production channel qq — > H ++ H and decay if ±=t — > 
liruf. The mass limits M H ±± > 110 150 GeV 0-Ei] were derived, with the strongest 
limits being for I = e, fx 15l-ll7j]. The branching ratios (BRs) for H ±zL 



depend on 

/; , , and are predicted in the HTM in terms of the parameters of the neutrino mass matrix 



14j, |19|, |20J. Detailed quantitative studies of BR(if 



(and BR(iP 



v)) in 



the HTM have been performed in [21H26I] with particular emphasis given to their sensitivity 
to the Majorana phases and the absolute neutrino mass i.e. parameters which cannot be 
probed in neutrino oscillation experiments. Simulations of the detection prospects of if ±=l= at 



the LHC with y/s = 14 TeV previously focussed on qq -> 7*, Z* H ++ H— only [27|, |28 
but recent studies now include the mechanism q'q — > H ±± H T [24], 25, 29]. The first search 
for H ±±: at the LHC with = 7 TeV [30| has recently been performed for both production 
mechanisms q'q — > H ++ H and qq' — > if ±± if =F , for the decay channels H ±i: — > if if and 
H ± — > ifu (where i,j = e, fx, r). 

In the HTM there are two electrically neutral mass eigenstates which are CP-even scalars. 
These are denoted by H\ and H2, with Mh 1 < Mh 2 - One of the eigenstates is dominantly 
composed of the isospin doublet field (and plays the role of the SM Higgs boson) while the 
other eigenstate is dominantly composed of the real part of the neutral triplet field. The 
mixing angle is small because it depends on the small ratio v&fv < 0.03 (where v = 246 GeV, 
the vev of the doublet field) . In phenomenological studies of the HTM it is usually assumed 
that the lighter eigenstate Hi is the one which is dominantly composed of the isospin doublet 
field. Therefore the phenomenology of Hi is more or less identical to that of the SM-Higgs 
boson. The converse case of the heavier eigenstate H 2 being the one which is dominantly 
composed of the isospin doublet field is possible in the HTM, and has been mentioned in 



3lH33j. However, no detailed study of the phenomenology of if 2 in such a scenario has 
been carried out. Importantly, if Mh 2 > [33| (where is one of the dominantly 

triplet eigenstates Hi, A , H^, if ±± ) then new decay channels for H2 would be possible. 
This would give rise to a phenomenology of H2 which differs somewhat from that of the SM 
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Higgs boson. In this work we focus on the case of Mh 2 > 2M H ±±, because a new production 
mechanism for H ±zL would be possible, namely gluon-gluon fusion gg — > H 2 followed by 
decay H 2 — > H ++ H . The case of Mh 2 > 2M H ±± necessarily requires Mh 2 ^200 GeV, a 
mass region which is now being probed for the first time by the LHC for the decay channels 



to SM particles, H 2 -> WW and H 2 -> ZZ |3J, |35 



Our work is organised as follows. In section II we describe the theoretical struc- 
ture of the HTM. In section III the parameter space for Mh 2 > 2M^ (where <fi is one 
of H\, A , H^, if ±± ) is described. In section IV the formulae for the decay widths of 
H 2 —> H ++ H , H + H~ , A A , H1H1 are presented. Section V contains a numerical analysis 
of the magnitude of the branching ratios of the above channels, as well as a quantitative study 
of the cross section for pair production of if ±± via gg — > H 2 , with decay H 2 — > H ++ H . 
Conclusions are given in section VI. 



II. THE HIGGS TRIPLET MODEL 

In the HTM @, @, a I = 1, Y = 2 complex SU(2) L isospin triplet of scalar fields 
is added to the SM Lagrangian. Such a model can provide a Majorana mass for the ob- 
served neutrinos without the introduction of a right-handed neutrino via the gauge invariant 
Yukawa interaction: 

C = h i:j tpf L CiT 2 Ai' jL + h.c (1) 

Here h^{i,j = e,fi,r) is a complex and symmetric coupling, C is the Dirac charge conjuga- 
tion operator, r 2 is a Pauli matrix, ipn = {i>i,ti) T L is a left-handed lepton doublet, and A is 
a 2 x 2 representation of the Y = 2 complex triplet fields: 

, 5+/V2 S ++ 
A= (2) 

A non-zero triplet vacuum expectation value (A ) gives rise to the following mass matrix 
for neutrinos: 

rriij = 2h ij (A } = V2hijV A (3) 

The necessary non-zero va arises from the minimisation of the most general SU(2) ® U(1)y 
invariant Higgs potential, which is written as follows 19|, |20j (with $ = (</> + , (f)°) T ): 



V(H, A) = -m%H^H + X(H^H) 2 + M| TrA f A + {^H T ia 2 A ] H + h.c.) 

+ Ai (if t i/)TrA t A + A 2 (TrA^) 2 + A 3 Tr(A t A) 2 + A 4 # f AA^H (4) 

Here m 2 H < in order to ensure (0°) = v/\^2 which spontaneously breaks SU(2) (£)U(1)y 
to U(1)q, and M\ (> 0) is the mass term for the triplet scalars. In the model of Gelmini- 



Roncadelli 36( the term fi(& T ir 2 A^&) is absent, which leads to spontaneous violation of 
lepton number for M A < 0. The resulting Higgs spectrum contains a massless triplet 
scalar (a Majoron, J) and another light scalar (H°). Pair production via e + e~ — > H°J 
would give a large contribution to the invisible width of the Z boson and this model was 
excluded at the CERN Large Electron Positron Collider (LEP). The inclusion of the term 
fi(§ T iT 2 A^) explicitly breaks lepton number when A is assigned L = 2, and eliminates 
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the Majoron [1, @, 0] . Thus the scalar potential in eq. (J4]) together with the triplet Yukawa 
interaction of eq. ([1]) lead to a phenomenologically viable model of neutrino mass generation. 
For small Va/v, the expression for va resulting from the minimisation of V is: 

va — —p= (5) 

V2(M| + ^(A 1 + A4)/2) 

For large Ma compared to v one has va — pv 2 /\/2M A , which is sometimes referred to as 
the "Type II seesaw mechanism" and would naturally lead to a small va- Recently there has 
been much interest in the scenario of light triplet scalars (M « v), (especially the distinctive 
doubly charged scalar, i? ±± ), within the discovery reach of the LHC, for which eq. (J5J) leads 
to Va ~ jU. In extensions of the HTM the term p,(<& T iT2A^<&) may arise in various ways: i) it 



can be generated at tree level via the vev of a Higgs singlet field [37( ; ii) it can arise at higher 



orders in perturbation theory [20j; iii) it can originate in the context of extra dimensions 



191 ] ; iv) it can arise in models with an additional heavy scalar triplet |38|. 

Some phenomenological studies focus on a simplified scalar potential (e.g. Ref. (24j ) 
in which the quartic couplings A,; (where i = 1,2,3,4) involving the triplet field A are 
neglected. The resulting scalar potential then depends on four parameters (—m 2 H , A, /i, 
Ma), but only three parameters are independent because the VEV for the doublet field 
{v = 246 GeV) is fixed by the mass of W^. The three independent parameters are usually 
chosen as A, Va, Ma or A, v a, p> (see eq. (J5J)). The inclusion of A* generates additional trilinear 
and quartic couplings among the scalar mass eigenstates. The terms with Ai and A4, which 
involve both triplet and doublet fields, are of particular interest because they can give a 
sizeable contribution to the masses of the scalar eigenstates (see below). A detailed study 
of the theoretical constraints on the scalar potential (e.g. vacuum stability, unitarity and 



perturbativity) has been carried out in [33 



An upper limit on Va can be obtained from considering its effect on the parameter 
p(= M^//M| cos 2 9yy). In the SM p = 1 at tree-level, while in the HTM one has (where 
x = va/v): 

1 + 2r 2 

The measurement p ~ 1 leads to the bound va/v <^0.03, or va ^8 GeV. Production 
mechanisms which depend on va (i.e. pp —> W^* —> W T H ±:iz and fusion via VF ± *W ± * — > 



if ±=t 12|, |39|) are not competitive with the processes qq —> H ++ H and qq' —> if ±± if =F at 
the energies of the Fermilab Tevatron, but such mechanisms could be the dominant source 
of F ±± at the LHC if v A = O (1 GeV) and M H ±± > 500 GeV. At the 1-loop level, v A must 
be renormalised and explicit analyses lead to bounds on its magnitude similar to the above 



bound from the tree-level analysis, e.g. see |40|, |41 



The scalar eigenstates in the HTM are as follows: i) the charged scalars iJ ±=t and H ± ; ii) 
the CP-even neutral scalars Hi and if 2; iii) a CP-odd neutral scalar A . The doubly charged 
if ±± is entirely composed of the triplet scalar field A ±=t , while the remaining eigenstates are 
in general mixtures of the doublet and triplet fields. However, such mixing is proportional 
to the triplet vev, and hence small even if va assumes its largest value of a few GeV. 1 



1 A large mixing angle is possible in the CP-even sector provided that Mh 1 ~ Mjj 2 
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III. SCENARIO OF A HEAVY SM-LIKE HIGGS BOSON (H 2 ) IN THE HTM 



In the HTM there are two CP-even mass eigenstates, which we denote by Hi and H 2 
(where both M H2 > M Hl and M Hl > M H2 are possible). Their compositions in terms of the 
original fields of the Lagrangian are as follows: 

H 2 = cos# h° + sm6 A H 1 = -sin# h° + cosfl A° (7) 

Here h° is the real part of the electrically neutral doublet field <p°, and A is the real part 
of the electrically neutral triplet field 5°. The mixing angle 9q is very small, 2 being of order 
0.03 at most (sin#o ~ v&fv). Hence H 2 is essentially composed of the doublet field h°, with 
couplings to the fermions and gauge bosons which are almost identical to those of the SM 
Higgs boson, while Hi is mainly composed of the triplet field A . 

The explicit expression for the 2x2 CP-even scalar mass matrix for the scalar potential 



in eq. Ml) is given in several works e.g. |3lH33lj. Neglecting the small off-diagonal elements 



in this mass matrix, the approximate expressions for the squared masses of Hi and H 2 are 
as follows: 

M 2 H2 = 2Xv 2 (8) 
M 2 Hl = Ml + + + 3(A 2 + A 3 K (9) 

The squared mass of the (dominantly triplet) CP-odd A is given by: 

M 2 A o = M 2 A + + ^)v 2 + (A 2 + X 3 )vl (10) 

The squared mass of the (dominantly triplet) H ± is given by: 

M 2 H± = Ml + + + (A 2 + V2X 3 )vl (11) 

Finally, the squared mass of the (purely triplet) doubly-charged scalar (if ±± = 5 ±± ) is given 
by: 

M 2 H±± = Ml + yi; 2 + X 2 v 2 A (12) 

One can see that the squared mass of the (dominantly doublet) H 2 is simply given by 2Xv 2 , 
as in the SM. In the expressions for the masses of M^ , Mg , M^ ± and M^ ±± there is a 
common term Ml + ^v 2 . It is evident that the mass scales for H 2 and the dominantly 
triplet scalars (A , Hi, H ± , if ±=t ) are unrelated, the former being set by 2Xv 2 and the latter 
by Ml + y v 2 ■ Neglecting the terms which are proportional to the small parameter v a, one 
can see that there are only two possible mass hierarchies for the triplet scalars, with the 
magnitude of the mass splitting being controlled by A4 (and M^o = MjjA: 

M A o,M Hl < M H ± < M H ±± for A 4 < (13) 
M H ±± < M H ± < M A o,M Hl for A 4 > (14) 



2 



The mixing angle can be maximal in the region of degeneracy Mh 2 ~ Mh 1 , but it quickly becomes small 
(~ va/v) with increasing mass splitting \Mh 2 — Afffi 
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In studies of the HTM it is sometimes assumed that M\ ^> 2\v 2 i.e. 
Mjj±±, Mjj±, Mhx, Mao S> Mff 2 . The motivation for this scenario is to have a "seesaw 
type" explanation for the smallness of va in eq. (JSJ). However, for Ma much larger than the 
TeV scale there would be no hope of observing the triplet scalars at the LHC In recent years 
there has been much interest in the study of the HTM as a TeV scale model of neutrino 
mass generation [ToT i.e. not invoking a large mass scale for M&. In these studies 



it is assumed (either explicitly or implicitly) that M\ > 2Xv 2 , with Ma< 1 TeV. 

The converse case where M A + ^-v 2 < 2\v 2 is rarely considered. In^ 3lM33| the possibility 



of Mh 2 > Mh±± , Mh± , Mn x , M40 has been mentioned, and in [33| the case of Mn 2 > 
2M^ (where <fi is one of H ±± ,H ± ,Hi,A°) is explicitly discussed. However, in these works 
there is no study of the phenomenology of H 2 at hadron colliders for the case of M# 2 > 
2M ( j ) , and how its experimental signature might differ from that of the SM Higgs Boson. 
Importantly, if M^ 2 > 2M ( f ) then new decay channels for H 2 become possible, 3 namely 
H 2 -> H ++ H—, H + H~, HxHx, A A . 

In this case the phenomenology of H 2 in the HTM could be different to that of the SM 
Higgs boson, because the new decay channels (if open kinematically) would compete with 
the usual decays of H 2 to SM particles (i.e. H 2 — > WW, ZZ, it). Of particular interest is the 
decay H 2 — > H ++ H , for which the condition M^ 2 > 2Mjj±± is necessary. If its branching 
ratio were sizeable then the production of H 2 via gluon-gluon fusion gg — > H 2 followed by 
the decay H 2 — > H ++ H would be an additional way to produce a pair of i? ±=t at hadron 
colliders. 

We note that the condition Mh 2 > 2M#±± necessarily requires Mh 2 <^ 200 GeV in order 
to respect the current lower bounds on M#±± from direct searches. At first sight, such a 
heavy SM-like H 2 would appear to be in conflict with experimental data, since it is well 
known that the Higgs boson in the SM is expected to be lighter than 200 GeV in order 
not to give an unacceptably large contribution to electroweak precision observables. In the 
context of the SM the case of Mh 2 ^ 200 GeV is quite strongly disfavoured, although this 



fact has not dissuaded direct searches in this mass region at the LHC [3J, [35j. However, the 
bound Mh 2 200 GeV cannot strictly be applied to the HTM, due to the additional scalar 
particles and the different renormalisation procedure, the latter being necessary because 
of the presence of the triplet vev (va). Dedicated analyses in models with scalar triplets 
have shown that a heavy (up to 1 TeV) SM-like Higgs boson can be made consistent with 



electroweak precision measurements [40, |41|. These studies are for a model with a real 



Y = scalar triplet, which has no doubly charged scalar and gives p > 1 at tree level, in 



contrast to the HTM which has p < 1 at tree level, (eq. ([6])). One can see in [4JJ that the 
condition Mh 2 > 2M^ (where <fi is one of the Y = triplet scalars) can be be accommodated. 
Although there is no explicit study in the HTM, we expect this result to also hold due to 
its greater number of free parameters (i.e. particle masses). In our numerical analysis we 
will treat M# 2 < 700 GeV and M^ 2 > 2M#±± as permissible parameter space in the HTM. 

From a phenomenological point of view, a heavy (>> 200 GeV) SM-like Higgs boson is 
attractive because it would be discovered more quickly at the LHC than a light SM-like Higgs 
boson with mass < 140 GeV. The region of 200 GeV < M# 2 < 500 GeV, for which the decays 
H 2 — > ZZ and H 2 — > WW are dominant in the SM, is a mass range where the LHC has 
sensitivity to cross sections which are much smaller than that of the SM Higgs boson. The 



The scenario of a heavy SM-like Higgs boson decaying to singly charged scalars, h° —> H + H , has been 



discussed in the Two Higgs Doublet Model 42 
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first searches at the LHC for a SM Higgs with Mh 2 > 200 GeV have already been carried 
out. The ATLAS collaboration (with 36 pb _1 of integrated luminosity) has searched for 
H 2 -> ZZ with the decay modes ZZ ->■ g+trw, ZZ -> and ZZ -> as well 

as if 2 — >• WW with the decay mode W W — > iisq'q [34| . The CMS collboration has searched 
for H 2 — >■ WW with the decay mode WW — >■ £z/£z/ [35(. Production of Hi is assumed to be 
via gluon-gluon fusion, gg — >■ i^2, and cross sections which are an order of magnitude above 
the prediction of the SM are currently being excluded at 95% c.l. By the end of the y's = 7 
TeV run (in which a few fb _1 of integrated luminosity will be accumulated), the sensitivity 
in these channels will be sufficient to exclude or provide evidence for the SM Higgs boson 
at a high confidence level in the region 200 GeV < Mjj 2 < 500 GeV. If the branching ratios 
of H2 —> H ++ H ,H + H~,HiHi,A°A° were sizeable then discovery of H2 in the channels 
H 2 — > ZZ and H 2 — > WW would require more integrated luminosity. 

IV. THE DECAYS H 2 -> # WITH = # ±=t , H ± , Hi, A 

There are four decay channels of H2 to pairs of scalars in the HTM: H 2 —> H ++ H , 
H 2 -»■ H+H~, H 2 A A and H 2 H X H X . If M H2 > 2M (where = H ±:t , i/ ± , A ) 
one can treat this as a two-body decay to a pair of on-shell <p. If M^ 2 < 2M^ we consider 
the partial width to be zero. Between one and four of the decays H 2 — > <p<f) can be open 
kinematically, depending on the mass splitting among <fi (which is controlled by A4 in eq. (fT3"j) 
and eq. 014p . For A4 > the lightest of the triplet scalars is if ±=t . In this scenario, once 
values of M# 2 and M#±± are chosen such that Mh 2 > 2M#±± there will be a value of A4 
above which only H 2 — > H ++ H is kinematically open. This will be the scenario where 
BR(# 2 -> H ++ H—) is maximal. 

The Feynman rules for the scalar trilinear couplings which mediate the decays are as 
follows (omitting a factor of —i): 



Here we consider H 2 to be entirely composed of the isospin doublet scalar field, which is true 
to a very good approximation. One can see that C^ 2 h++h — is controlled only by Ai, while 
the other trilinear couplings depend on both Ai and A4. If Ai and A4 are sizeable, then the 
branching ratios for H 2 ->■ H ++ H , H 2 ->■ H + H~, H 2 ->■ H\H X , and H 2 -)> A A could be 
non- negligible. 

One can use a generic formula for the decay rate for the four channels: 



Here 5h = 2 for <p = i? ±=t , (i.e. non-identical particles in the final state) and 5 = 1 for 
(f) = Hi, A (i.e. identical particles in the final state) 

It is clear that the two crucial parameters for a large BR(if 2 — * H ++ H ) are Ai 
(which determines the strength of the coupling \Ch 2 h++h — I in eq. (ITS"]) and M H ±± (which 
determines the suppression from phase space). In our numerical analysis we shall take 



Cr 2 h++h— = Aiu 

„ / , A4 . 

C H2 h+h- = (Ai + —)v 

Ch 2 h 1 h 1 ,Ch 2 a ( >a ( > — (Ai + A 4 )u 



(15) 
(16) 
(17) 




(18) 
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Mh±± as an input parameter. As can be seen from eq. (I12p . the dominant contribution 
to Mjj ±± is from the two terms M\ + \\v 2 /2. Therefore, by taking Ai and M#±± as in- 
put parameters the value of M\ is determined. We will be focussing on the parameter 
space of 90 GeV < M H ±± < 300 GeV and < Ai < 4, and consequently M\ < when 
Mfj ±± < X 1 v 2 /2. In the scenario of M\ < the positive mass of M#±± is obtained from the 
term Xiv 2 /2. Alternatively, one could consider Ai < and M\ > 0. The crucial point here 
is that the parameters M\ and Ai should have opposite signs if one wishes to have large |Ai| 
(in order to enhance \Ch 2 h++h — I) together with a fairly light if ±=l= . 

We now summarise the current lower limits on M#±± from direct searches. There have 
been searches for the decay channels H ±i: — > if if for i,j = e, u,r (these are the dominant 
decay channels for v A £0.1 MeV) at LEP 0, Tevatron ^Bum and the LHC (CMS 
Collaboration). The strongest mass limits for the decays H ±:L — > ee,efi,fifi are from the 
LHC search, which obtained M#±± > 144, 154, 156 GeV respectively, assuming BR= 100% 
in a given channel. Separate searches for three and four leptons (which have significantly 
different backgrounds) were performed. These limits are weakened considerably for the case 
of BR< 100% because the event number for the signal is proportional to the square of 
the branching ratio (BR 2 ). In [3(| both production mechanisms qq — > 7*,Z* — > H ++ H~ 
and q'q — > W — > H ±:iz H T were considered (with the assumption M#±± = M#±), which 
increases the sensitivity in the three-lepton channel. For the decays involving one r, namely 
J ff ±± er,/xr, the limit M H ±± > 106 GeV was derived in both channels at the LHC 



30|, with stronger limits from the Tevatron (M#±± > 112 GeV and 114 GeV respectively) 
obtained in 18j. The only search for if ±± — > r ± r ± at a hadron collider is the LHC search 



m 



30J, which derived the limit M H ±± ^80 GeV. 
We will take M H ±± = 90 GeV as our lowest value for the mass of -ff ±=t , and this is 
allowed for certain choices of branching ratios of i? ±± . As explained above, the limits on 
Mh±± from hadron colliders are weakest for those channels which involve r. In contrast, the 
limit from the LEP searches of M#±± ;> 100 GeV applies to all the decays if ±± — > if if with 
i, j = e, yU, r. The search strategy at LEP requires four leptons and so the event number for 
the signal is proportional to BR 2 . The scenario of M#±± = 90 GeV is compatible with the 
all the above searches provided that the decays involving r are dominant e.g. choices like 
BR(if ±=t — > eT,fiT,rr) of around 33%. It is not necessary to have BR(iJ ±=t — > ee,efi,fifi) 
totally absent for M H ±± = 90 GeV, and BRs of the order of 10% for these channels can 
be accommodated because the event number is proportional to BR 2 , and for BR = 10% 
this is a large suppression factor. We note that the sum of BR(i7 ±=b — > ee, e/i, /i/i) cannot 
be taken arbitrarily small in the HTM because the Yukawa couplings are related to the 
neutrino mass matrix via eq. ([3]). The allowed values of BR(i/ ±± ->■ if if) in the HTM 
have been studied in detail in |2l|-|26|. and in 25| it can be seen explicitly that the sum of 
BR(if ±=t — > ee,efi,fifi) must be greater than around 5%. 

Very recently the searches for i? ±=t — > ee, ejj^fifi by the CDF collaboration in [l5[ (which 



used 0.24 fb" 1 ) were updated using 6.1 fb" 1 |44j. Mass limits of M H ±± > 225, 210, 245 GeV 



were obtained, again assuming BR= 100%. In these searches the event number for the signal 
is linear in BR, and for BR~ 3%(15%) the limit M H ±± > 245 GeV for /i/i would 

weaken to Mn±± > 100 GeV (150 GeV). Note that these mass limits in 44j only assume 
production of if ±± from qq — > 7*, Z* — > H ++ H . The inclusion of q'q — > W — > H ±:L H T 
would allow larger values of M H ±± to be probed. Finally, if the decay H ±ziz — > W^W^ 1 is 
dominant (which is the case for va ^0.1 MeV) then M H ±± = 90 GeV is permitted because 
there have been no direct searches for this channel in the context of models with H ±:k . We 
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will respect the all the above mass limits in our numerical analysis, the most stringent ones 
being for the channels if ±=t — > e ± e ± , e fi and fi ± fi ± . 

We will only consider the scenario of Mh 2 > 200 GeV for which the decay channels 
H 2 — > WW and H 2 — > ZZ can be treated as two-body decays. The expressions for their 
partial decay widths are as follows: 

T(H 2 -> WW) = ^—^{l ~ + 124.) (1 - 4k w ) 1/2 <V (19) 

OZ7T 

with k w = M^/Mjj 2 and 5w — 2. 

T{H 2 ZZ) = F H2 (l - Ak z + 12k|)(1 - AK Z f/ 2 5 z (20) 
6Aix 

with k z = M z /Mjf and 5 Z = 1. If M H2 > 2m t then the decay channel H 2 — > it is open: 

T(H it) = M H fi (21) 

where /3 t = (1 — Am 2 t / M^) 1 / 2 . All other decays of if 2 to SM particles (e.g. H 2 ->■ 66, r+r") 
have negligibly small partial widths for M// 2 > 200 GeV. Note that other decay channels 
such as H 2 — > H ± W, H 2 — > H^W and H 2 — > A°Z are suppressed by the small mixing 
between the doublet and triplet fields, and so can be neglected. 



V. NUMERICAL ANALYSIS 

We now study the magnitude of the branching ratios of the decays channels H 2 — > <fxp 
for (f) = if ±=t , H^, Hi, A . The four important parameters are Mjj±±, Mh 2 , X\ and A4. 
The other parameters in the scalar potential are fixed as A2 = A3 = 0.5 and Va = 10 -2 
MeV, the latter choice ensuring that the decays if ±=t — > £ ± £ ± are dominant. These latter 
three parameters appear in the expressions for the masses of the triplet scalars in eq. (Q to 
eq. (I12p but their effect is essentially negligible, even for the case of v a = 1 GeV. We will 
present results for M H ±± = 90 GeV, 150 GeV, 200 GeV and 300 GeV. As explained in the 
previous section, the choice of M H ±± = 90 GeV requires small BRs (< 3%) for the decay 
channels i? ±=t — > e ± e ± ,e ± fi ± and /i ± / u ± (and consequently large BRs to channels involving 
r) in order to respect the limits from the direct searches for if ±± . Larger values (^> 3%) of 
BR(H ±zt — > e ± e ± ,e ± fi ± , ^fi^ are permitted as M H ±± increases. 

Fig. <^l) shows the branching ratios of H 2 as a function of M# 2 . We take M#±± = 90 
GeV, A = 1, and A 4 = 0.8 (the latter choice gives M H ± = 142 GeV and M A ^ H o = 179 GeV). 
The magnitude of BK(H 2 — > H ++ H ) can reach 65% for Mh 2 = 200 GeV, and stays as the 
dominant channel until Mh 2 ~ 260 GeV, at which H 2 — > WW becomes dominant. BR(i^2 - > 
H++H-) falls below BR(# 2 ZZ) at around M H2 = 320 GeV. This dependence on M H2 
can be explained by the fact that the partial widths of H 2 — > WW, ZZ are proportional to 
Mfj 2 , and so ultimately these channels will dominate for larger M# 2 . The other decays of 
H 2 to two triplet scalars also can have sizeable branching ratios, with BK(H 2 — > H + H~) 
reaching 20% at most, and exceeds B>R(H 2 — > H ++ H ) for Mh 2 ^315 GeV. The branching 
ratios of H 2 — > A A and H 2 — > H\H\ are equal; they are plotted individually and their sum 
peaks at ~ 10%. In fig. ((lb) we show contours of BR(ii2 — > H ++ H ) in the plane [M# 2 , AjJ. 
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As expected, BR(H 2 — > H ++ H ) takes its largest values for large Ai and light Mh 2 , with 
BK(H 2 — > H ++ H ) > 90% being possible. We note that such a scenario would render the 
searches for H 2 — > WW, ZZ ineffective until a very large amount of integrated luminosity is 
obtained. In fig. ([Tb) we show contours of BR(if 2 — > H ++ H ) in the plane [A 4 , Ai], fixing 
Mh 2 = 300 GeV. For A4 ;> 1 only the decay H 2 — > H ++ H is open kinematically and so 
the contours are horizontal. For A4 = all the triplet scalars are degenerate and thus all 
four decay channels are open. Figs. I2"f3l and H] are analogies of fig. [U but with M H ±± = 150 
GeV, 200 GeV and 300 GeV respectively. In all figures we take Mh 2 > 2M H ±± . One can see 
similar qualitative behaviour, but since the lowest value of Mh 2 is larger in figs. 121131 and H] 
than in fig. [H the maximum values of BR(iJ 2 — > H ++ H ) are less than in fig.l. However, 
in fig. 2b, fig. 3b and fig. 4b it can be seen that BR(# 2 -> H ++ H—) > 50%, 25%, 5% 
respectively is possible for Ai <^3. 

It clear that BR(if 2 — > H ++ H ) can be sizeable, and we will now quantify the magni- 
tude of the pair production of H ±i: which originates from production and decay of H 2 . 
At hadron colliders H 2 is dominantly created via gluon-gluon fusion, gg — > H 2 . For 
Mh 2 = 2M H ±± the cross section of gg — > H 2 at the LHC is significantly larger than the 
cross section for the direct production mechanisms of H ±i: (i.e. qq — > 7*,Z* — > H ++ H~ 
and q'q — > W — > H ±:k H^). However, the same is not true at the Tevatron, and 
a(gg -> H 2 ) <,a(qq 7*, Z* H ++ H") for M H2 > 2M H ±±. 

We introduce the ratio R, defined by: 

R = o(99 H 2 ) x BR(ff 2 -» 

<r(qq->'f,Z* -> H++H-) 1 ' 



The denominator in eq. ( 1221) is the conventional mechanism for production of H ++ H , 
which is assumed in the ongoing searches for if ±=t . The numerator is a novel mechanism 
which contributes when BR(/7 2 - > H ++ H ) 7^ 0. We will now study the magnitude of the 
ratio R at the LHC (with yfs = 7 TeV and 14 TeV) and at the Tevatron. In Fig. [5] we plot 
R as a function of M H2 at the LHC with y/s = 14 TeV, for M H ±± = 90 GeV, 150 GeV, 
200 GeV and 300 GeV. The factorisation scale and normalisation scale are both taken to 
be Mh 2 for gg — > H 2 , while for qq — > 7*, Z* — > H ++ H both scales are taken to be the 
partonic centre-of-mass energy. We use CTEQ6L1 parton distribution functions [45[ with 
the leading-order partonic cross section for gg — > H 2 [46|. We do not apply QCD K factors 
which, would increase the value of R because the ratio of the K factors for a(gg — > H 2 ) 



47j and a(qq -»■ 7*, Z* — >■ H ++ H ) [10( is about 1.4 in the region of interest of M H2 and 
M H ±±. 

In fig. we take Mu±± = 90 GeV, which fixes the value of cr(qq — > H ++ H ), and 
A 4 = 0.8. We take A x = 1 and 4. If \ x = 1 one can see that R = 4.7 for M H2 = 200 GeV, 
and R > 1 for M H2 < 290 GeV. If Ai = 4, one finds that R = 7.0 for M H2 = 200 GeV, and 
R > 1 for M H2 < 420 GeV. The noticeable drop in the value of R for M Hi ~ 280 GeV is 
due to the opening of the decay channel H 2 H + H~ (see fig. (Jit))- Both cr(gg — > H 2 ) and 
BR(if 2 — > H ++ H ) are decreasing functions of M# 2 , which explains the overall dependence 
of R on Mh 2 - Note that R does not fall so sharply with Mjj 2 i n the region 320 GeV < 
M H2 < 380 GeV, because cr(gg —> H 2 ) increases in magnitude up to a local maximum at 
M H2 = 2m t , before decreasing again. In fig. ([5b) we take M H ±± = 150 GeV, and R ~ 16 
for Mh 2 = 2m t and Ai = 4. Larger values of R are attainable because the magnitude 
of <j{qq 7*,Z* — > H ++ H ) (i.e. the denominator eq. fT22|) ) diminishes considerably 
when going from Mh±± = 90 GeV to Mh±± = 150 GeV, while the corresponding decrease 
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in a(gg — > if 2) for larger Mh 2 is relatively less. In fig. (j5b) (for M#±± = 200 GeV) the 
maximum value is i? ~ 19, and in fig. (j5H) the maximum value is R ~ 4. It is evident that 
there is a sizeable parameter space for R > 1, and thus gg — > H 2 could give a significant (or 
even dominant) contribution to the pair production of H ±=l= at the LHC. We also note that 
the decay H 2 — > H + H~ (which can have a large BR in fig. ([Th.) — > fig. (|3J0) can l ea d to 
additional production of H ++ H because the branching ratio of the decay H ± — > H^W* 
can be large in a sizeable parameter space of [ua ; M h ± — M H ±±], as shown in (48[. In fact, 
in fig .(UK) fig. (H» the mass splitting (M H ± - M H ±±) is between 20 GeV and 52 GeV, 
and with our chosen value of v& = 10~ 2 MeV one would have BR(if ± — > H ±:k W*) > 99%. 

In fig. ([6]) we plot the analogies of fig. ([5]) for the LHC with y/s = 7 TeV. One sees 
a similar qualitative behaviour, with lower maximum values of R. In fig. ([7]) we plot the 
corresponding results for the Tevatron, for M#±± = 90 GeV and M#±± = 150 GeV. Since 
a(gg -> H 2 ) ^<r(qq 7*, Z* H ++ H ) for M H2 > 2M H ±± at the Tevatron, the maxi- 
mum value of R ~ 0.4 (for Ai = 4) is much smaller than at the LHC and is comparable to 
the QCD K factor for a(qq -> 7*, Z* -> H ++ H~~) [h}. 

Finally, we quantify the number of H ++ H events for a given integrated luminosity £ 
at the LHC. We introduce the parameter Nh±±, which is defined as follows: 

N H ±± =exCx [a{qq 7*, Z* -> H ++ H—) + a(gg -> Jf 2 ) x BR(# 2 -> H ++ H~~)] (23) 

The efficiency e is the fraction of H ++ H events which remain after all acceptance/selection 
cuts are imposed to reduce the background from the SM. The value of e depends on which 
decay channel if ±± — > if if is being considered. From the LHC simulation in j28j for the 
decay if ±=t — > with y/s = 14 TeV, one can derive an approximate value of e w = 0.73 

for M H ±± = 600 GeV and e MM = 0.64 for M H ±± = 300 GeV. As expected, the efficiency 
is greater for larger M^±±, since the leptons originating from H ±±: are more energetic. 
Extrapolating these values to the region of M H ±± < 300 GeV (the mass region on which 
we will focus) one roughly obtains 0.5 < e MAt < 0.6. The efficiencies for the decay channels 



if ±=t — > e ± e ± and H ±±L — > e ± /j. ± are expected to be similar in magnitude to e MM (see |30|). 
The efficiencies for the decays of H ±i: to final states involving a r lepton are much lower 
e.g. in 3(| one can derive e Mr ~ 0.02 for the channel if ±=l= — > /i ± r ± , with even lower values 
for the channel if ±=l= — > t ± t ± . We will show results for the decay mode H ±i: — > \i \i , for 
M H ±± = 90 GeV, 150 GeV, 200 GeV and 300 GeV. The strongest lower bounds on M H ±± in 
this channel (assuming a branching ratio of 100%) are M#±± > 156 GeV from the LHC in 
r 30( | and Mh±± > 245 GeV from the Tevatron in 44], both limits being preliminary results. 



For the case of BR(if ±=t — > fi^fJ^) < 100%, one can derive from 144] the approximate 
limits M H ±± ^100 GeV, ^150 GeV and ^200 GeV for BR £3%, £15% and £40%, 
respectively. We do not include these values of BR(if ±=t — > yU ± /i ± ) when showing results for 
Nh±±. In future searches which require three or four leptons (as done in the LHC search 
in j30|) the event number Nh±± in eq. (1231) needs to be scaled by a multiplicative factor of 
[BR(if ±=t — )• if if)] 2 - As explained above, this factor is necessarily less than unity if one 
considers Mh±± < 245 GeV with decay H ±i: — > /j^/j^. 

In fig. © we show N H ±± in the plane [M H2 , Ai] for y/s = 14 TeV with £=30 fb -1 . We 
use e w = 0.64 for M H ±± = 300 GeV, and e MA1 = 0.50 for the other chosen values of M H ±± 
(90 GeV, 150 GeV, 200 GeV). The contribution to N H ±± from qq -> 7 *,Z* H ++ H~ 
alone does not depend on M H2 and Ai, and is roughly equal to 20500, 3270, 1130 and 280 for 
M H ±± = 90 GeV, 150 GeV, 200 GeV and 300 GeV respectively. In each panel in fig. © the 
contour with the lowest number of events corresponds to a value of Nu±± which is slightly 
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larger than the above values for Nh±± from qq — > 7*, Z* — > H ++ H alone. We emphasise 
that the displayed Nh±± for M#±± = 90 GeV, 150 GeV and 200 GeV need to be multiplied 
by the square of BR (for a future three or four lepton search) where BR ^3%, <> 15% and 
<^ 40% in order to comply with the mass limits in jZIJ . Fig. OH]) can also be applied to other 
decay channels such as H ±±L — > fi^r^ 1 after multiplying N H ±± by e^ T /e m ~ 1/30. Clearly, 
the contribution from a(gg — > H 2 ) x BR(if 2 — > H ++ H ) could significantly enhance the 
number of H ++ H events at the LHC, provided that M^ 2 > 2M H ±± and Ai is not very 
small. Since it is not expected that Mh 2 ^ 700 GeV (from considering constraints from 



perturbativity and unitarity e.g. see 33]), the enhancement from a(gg — > H 2 ) x ~BR(H 2 



H ++ H ) is limited to the region M#±± <; 350 GeV. However, its contribution would allow 
the possibility of probing smaller values of BR(if ±± — > if if) for a given value of M#±± 
(provided that M H2 > 2M H ±±). 

In fig. §S§ we show N H ±± for \/s = 7 TeV with C=2 fb -1 . We take a slightly lower 
efficiency of e MAt = 0.4 for M H ±± = 90 GeV, 150 GeV and 200 GeV, which is in rough 
agreement with the value for the channel H ±i: — > [i^y^ in the experimental search at 
y/s = 7 TeV in {30[. For M H ±± = 300 GeV we take e MA1 = 0.5. Again, the enhancement from 
a (99 ~^ H2) x BR(i^2 - > H ++ H ) can be sizeable, and could lead to a quicker discovery 
of a light H ±i: with BR(H ±i: if if) < 100%. We do not plot an analogous figure for 
the Tevatron since the maximum value of R is much smaller than at the LHC, as shown in 
fig- ©. 

Finally, we emphasise that the parameter space of Mh 2 > 2M H ±± will be probed by two 
distinct searches with the LHC data taken at y/s = 7 TeV: i) the search for H 2 —> WW, ZZ 
(with first results in [Hi EHf), and ii) the search for qq — > 7*, Z* — >■ H ++ H (with first 



results in [30|, and a simulation in |49[). The run at y/s = 7 TeV, with possibly up to 5 
fbr 1 of integrated luminosity, has the potential to exclude or provide evidence for a SM-like 
Higgs boson with 200 GeV < M^ 2 < 500 GeV at a high confidence level. Therefore the 
scenario of M^ 2 > 2M#±± and its possible impact on the direct searches for i7 ±=t should be 
clarified within the next two years. 



VI. CONCLUSIONS 

Doubly charged Higgs bosons (if ±=t ), which arise in the Higgs Triplet Model (HTM) 
of neutrino mass generation, are being searched for at the Tevatron and at the LHC. 
The ongoing searches assume the production mechanisms qq — > 7*, Z* — > H ++ H and 
q'q — > W* — > f/" ±± i/ =F , with the leptonic decay H ±±L — > ifif- We proposed an additional 
production mechanism for i? ±=t , which becomes possible if the SM-like Higgs boson in the 
HTM (H 2 ) is heavy enough to decay to a pair of on-shell if ±=t . We quantified the magni- 
tude of the branching ratio of H 2 —> H ++ H , and showed that it can be large (>> 10%) 
if a quartic coupling in the scalar potential is sizeable, Ai > 1. We performed a numerical 
study of the production rate for H 2 via gluon-gluon fusion, gg — > H 2 , followed by the decay 
H 2 — y H ++ H , and we showed that its cross section at the LHC can be greater than that of 
qq — > 7*, Z* — > H ++ H in a sizeable parameter space of [Ai, M H2 ] (with M H2 > 2M H ±±). 
In the optimal case (e.g. Ai = 4, M H2 ~ 420 GeV and M H ±± = 200 GeV) the ratio of 
cross sections can be as large as nineteen. We note that our analysis was carried out using 
the leading-order cross sections only, and the inclusion of QCD K factors would provide 
a further enhancement of 40% in the above ratio. Such an additional source of H ++ H~ 
would enable smaller values of the branching ratio of if ±=t — > if if to be probed at the LHC. 
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FIG. 1: Upper left panel a): branching ratios of H 2 as a function of Mjj 2 . Upper right panel b): 
Contours of BR(if2 — > H ++ H ) in the plane [M# 2 , Ai]. Lower panel c): Contours of BR(if2 — > 
H+ + H in the plane [A4, Ai]. In all figures Mjj±± = 90 GeV. In a) and b) A4 = 0.8, which gives 
M H ± = 142 GeV and M A ^ Hl = 179 GeV. In c) M H2 = 300 GeV. 



The case of Mh 2 > 2M H ±± necessarily requires Mh 2 ^ 200 GeV, and this mass region is now 
being probed for the first time at the LHC for the decay channels of H 2 to SM particles, 
H 2 ->■ WW and H 2 -> ZZ. The possibility of a SM-like Higgs boson in the HTM with 
Mh 2 > 2M H ±± and its potential impact on the direct searches for H ±ziz should be clarified 
within the y/s = 7 TeV run at the LHC. 
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FIG. 2: Upper left panel a): branching ratios of H 2 as a function of Mjj 2 . Upper right panel b): 
Contours of BR(if2 —> H ++ H ) in the plane [M# 2 , Ai]. Lower panel c): Contours of BR(i^2 — > 
H ++ H in the plane [A4, Ai]. In all figures M H ±± = 150 GeV. In a) and b) A4 = 0.8, which gives 
M H ± = 186 GeV and M A ^ Hl = 216 GeV. In c) M H2 = 400 GeV. 



Note Added 



After submission of this paper, the LHC searches for the SM Higgs boson were up- 
dated with C = 1.1 fb" 1 0j. For the region M H2 > 200 GeV, both of the CMS 
and ATLAS collaborations use the decay channel H 2 — >• ZZ, with subsequent decays 
ZZ £ + £-uu,ZZ -»■ e + e~qq and ZZ -> £+£-£+£-. CMS also search for H 2 ->■ VFVF 
with the decay mode H^VF — > &v£v, while ATLAS search for H% — >■ VFVF with the decay 
mode VFiy — > £vq'q. After combining the results from these four distinct channels, both 
collaborations exclude at 95% c.l the mass range 295 GeV < Mh 2 < 450 GeV. This does not 
preclude a sizeable value of R in the HTM, e.g. from fig. [6fc, one can see that 7 > R > 1 in 
the interval 450 GeV < M# 2 < 600 GeV, for Ai = 4 and y/s = 7 TeV (and not including the 
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FIG. 3: Upper left panel a): branching ratios of if 2 as a function of Mjj 2 . Upper right panel b): 
Contours of BR(i^2 — > H ++ H ) in the plane [M# 2 , Ai]. Lower panel c): Contours of BR(i^2 — > 
H ++ H in the plane [A4, Ai]. In all figures M H ±± = 200 GeV. In a) and b) A4 = 0.8, which gives 
M H ± = 228 GeV and M A ^ Hl = 253 GeV. In c) M Ha = 500 GeV. 



enhancement from the QCD K factor). 
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FIG. 7: The ratio R = a(gg H 2 ) x BR(H 2 -> H ++ H—)/a{qq -> H ++ H—) at the Tevatron 
(with ^/s = 1.96 TeV) as a function of M# 2 . The curves are for Ai = 4 and A4 = 0.8 . We take 
M H ±± = 90 GeV in panel (a), M H ±± = 150 GeV in panel (b). 
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FIG. 8: The number of H ++ H~ events (assuming BR(# ±:t -)• ^ ± /i ± ) = 100%) at the LHC with 
y/s = 14 TeV and £ = 30 fir 1 in the plane [M# 2 ,Ai]. We take M H ±± = 90 GeV in panel (a), 
M H ±± = 150 GeV in panel (b), M H ±± = 200 GeV in panel (c) and M H ±± = 300 GeV in panel 
(d). In all figures A4 = 0.8. The number of events for qq — > 7, Z — > H ++ H alone is 20500, 3270, 
1130 and 299 in (a),(b),(c) and (d), respectively. 
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FIG. 9: The number of H ++ H~~ events (assuming BR(H ±i: -)• ^ ± /i ± ) = 100%) at the LHC with 
= 7 TeV and C = 2 fbr 1 in the plane [M# 2 ,Ai]. We take M^±± = 90 GeV in panel (a), 
M H ±± = 150 GeV in panel (b), M H ±± = 200 GeV in panel (c) and M H ±± = 300 GeV in panel 
(d). In all figures A4 = 0.8. The number of events for qq — > 7, Z — > H ++ H alone is 390, 53, 16, 
3 in (a),(b),(c) and (d), respectively. 
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